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Abstract 

We analyse the extraction and the subsequent fall of a coherent atomic beam from a trapped 
Bose-Einstein condensate using a rf transition to a non-trapping state, at T=0 K. Our treatment 
fully takes gravity into account but neglects the details of the mean-field potential exerted on the 
free falling beam by the trapped atoms. We focus on the weak coupling regime, where analytical 
expressions for the output rate and the output mode of the "atom laser" , can be obtained. Com- 
parison with experimental data of Bloch et al. [Phys. Rev. Lett. 82, 3008 (1999)] without any 
adjustable parameter is satisfactory. 
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Bose-Einstein condensates (BEC) of dilute alkali vapors constitute a potential source 
of matter waves for atom interferometry, since it has been proven [Q] that they are inher- 
ently fist-order coherent. Various schemes for "atom lasers" have been used to extract a 
coherent matter wave out of a trapped BEC. Pulsed devices [0, ^, |^ were demonstrated by 
using an intense spin-flip radio frequency (rf) pulse, Raman transitions or gravity-induced 
tunneling from an optically trapped BEC. Later on, a quasi-continuous atom laser has been 
demonstrated by using a weak rf field that continuously couples atoms into a free falling 
state 0|. This "quasi-continuous" atom laser promises spectacular improvements in appli- 
cation of atom optics, for example in the performances of atom-interferometer-based inertial 
sensors 0. 

Gravity plays a crucial role in outcouplers with spin-flip rf transitions: it determines the 
direction of propagation of the extracted matter wave, its amplitude and its phase. How- 
ever, most of the theoretical studies do not take it into account (for an up-to date review, 
see Ref.p]). Although gravity has been included in numerical treatments ^ relevant for 
the pulsed atom laser of Mewes et al. H, to our knowledge, only the ID simulations of 



Refs.g, |T0| treat the quasi-continuous case in presence of gravity, and the results of ||TG 
compare only qualitatively to the experimental data of 0. 

In this Letter, we present a 3D analytical treatment fully taking gravity into account. The 
extraction of the atoms from the trapped BEC and their subsequent propagation under 
gravity are analytically treated in the weak coupling regime relevant to the quasi-continuous 
atom laser. We derive an expression for the atom laser wave function and a generalized rate 
equation for the trapped atoms that agrees quantitatively with the experimental results of 
1^, without any adjustable parameter. 

We consider a '^^Rb BEC in the F = 1 hyperfine level at T = K. The m =-1 state is 
confined in a harmonic magnetic potential Krap = + i^lz/^ + ^^^z^)- A rf magnetic 

field Brf = Srf cos (c(Ji.ft)ex can induce transitions to m =0 (non-trapping state) and m =+1 
(expelling state). The condensate three-component spinor wavefunction ^ = [%l)'^rn=-ifl,+i 



obeys a set of coupled non-linear Schrodinger equations [11 . Following Steck et al. [12 



we consider the "weak coupling limit". In this regime to be defined more precisely later, 
the coupling strength is low enough that the populations of the three Zeeman sublevels 
obey the following inequality: A^+i <^ A'^g <^ A^_i. In the rest of this paper, we therefore 
restrict ourselves to m =-1 and m =0, and set the total atomic density n{v) k,\ ilj_i[r,t) [^. 
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At this stage, the components ipm = V'm'e *™-'^rft obey, under the rotating wave approxima- 
tion, the following two coupled equations: 

= [Mrf + + ^v^o (1) 

with n^i = pV2M + Krap + f/ I ^-1 P and Hq = p^/2M-Mgz + U \ p. The strength 
of interactions is fixed by U = Anh^aN/M, a ~ 5 nm being the diffusion length that we 
suppose equal for all collision processes, and the initial number of trapped atoms. The 
rf outcoupler is described by the detuning from the bottom of the trap hS^-i = Vos — ^t^rf 
and the Rabi frequency h^lj-f = ii-qB^i /2\/2 (taking the Lande factor g-p=-l/2). The origin 
of the z axis is at the center of the condensate, displaced by gravity from the trap center by 
Zsag = g/^±- We have taken the zero of energy at 2; = in m =0, so that the level splitting 
at the bottom of the trap is Vos = /iB-Bo/2 + Mg"^ /2uj']_ (Bq is the bias field). 
The evolution of the m =-1 sublevel at T = K is described by a decomposition of 



the wave function into a condensed part and an orthogonal one that describes elementary 
excitations (quasiparticles). As H-i depends on time through | ip_i |, both the coefficients 
and the eigenmodes depend on time. However, we assume an adiabatic evolution of the 
trapped BEC [0, so that the uncondensed component remains negligible and: 

^_i(r,t) ^ «(t)0_i(r,t)e-'i'o(^(*')/'*+'5^f)'^*' (3) 

where a(t) = (A^_i(t)/A^)^/^. The time-dependent ground state of energy /i is, in 

the Thomas-Fermi (TF) approximation 0_i(r, t) = (yu/f/)-'^/^[l — — 5^]^^^, where 

f^=(x/xo)^ + {y/yoY and z=z/zq are such that f\ + z^ < 1. The BEC dimensions are 
respectively Xo=(2/i/Mu;^)^/^ and yQ=ZQ={2iJi / Mu"]}^/"^ . The condensate energy is given 
by ii={nuj /2)(lbaN^i/ af/^ (we have set ^={^^^1)^/^ and a={n/Mujf''^). The time- 
dependence of 0_i is contained in /x, xq, and zq, which decrease with A^_i(t). In the 
following, we will take typical values corresponding to the situation of (reported in [p^): 
N = 7 X 10^ atoms initially, uJx=2'k x 13 Hz and uj±=2tt x 140 Hz, which gives Xq ~ 55.7/im, 
zq ~ 5.3/im and fi/h ^ 2.2 kHz. 

For the m=0 state, the hamiltonian of equation (^ with =0 becomes in the TF ap- 
proximation [|l^: 7^0 =p^/2M — Mgz + /iMax[0, 1 — — 5^]. There is a clear hierarchy 
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in the energy scales associated with each term in TYq- First, in the TF regime, the kinetic 
energy in the trap is neghgible as compared to the mean field energy. Second, for the atom 
numbers involved in current experimental setups, fi/MgzQ is small, and the mean field term 
is only a perturbation as compared with gravity. The leading term is therefore the gravi- 
tational potential, that subsequently converts into kinetic energy along z. On the contrary, 
the transverse kinetic energy (along x and y) remains small (we will make use of this to 
compute the output rate). In addition, we neglect for simplicity the spatial variations of the 
mean field potential: we replace it by the mean-field contribution to the energy /i needed 
to remove one atom, i.e. = 4yu/7 @]. In summary, we consider in the following the 
approximated hamiltonian Hq ^ p^/2M — Mgz + 4/x/7. 

The eigenstates of Ho are factorized products of one dimensional wavefunctions along the 
three axis. In the horizontal plane (x,y), the eigenstates are plane waves with wavevectors 
kx,ky that we quantize with periodic boundaries conditions in a 2D box of size L. Conse- 
quently, the wavefunction is 0j}(a;, y) = L^'^ e^^^'^^^^^^'^ and the density of states p^y = L^/47r^. 
Along the vertical direction z, the normalizable solution of the ID Schrodinger equation in 
a gravitational field is 00^^^"*= ^■^i{~CEj, where Ai is the Airy function of the first 
kind and (e^ = ~ ^eJ/^- The classical turning point ze^ = —E^/Mg associated with the 
vertical energy labels the vertical solution; / = {h^ /2M'^gY/^ is a length scale, such that 
I <^ xo, 2/05 zo (for ®^_R6, / ^ 0.28/im). In order to normalize Ai and to work out the density 
of states, we restrict the wavefunction to the domain ] — oo, where the boundary at z=H 
can be arbitrarily far from the origin. In ] — oo, zeJ\, Ai falls off exponentially over a distance 
/, while it can be identified with its asymptotic form Ai{—s) ~ vr^/^s^"'^/'* cos(2s^/^/3 — vr/4) 
for s ^ 0. To leading order in iJ, by averaging the fast oscillating cos^ term to 1/2, we obtain 
the normalization factor A={tiH~^^'^ /ly/'^ . The longitudinal density of states follows from 
pz{E^) = h-^dV/dE, where V = J^^e, dzdp^e{pl/2M - Mgz- E^) (with 9 the step function) 
is the volume in phase space associated with energies lower than Ez . We restrict the z space 
to the domain defined above, and the space to p^ > 0. In this way, we do not count the 
component of Ai that propagates opposite to gravity|T7[]. We obtain pz{zEj = {l/2nl)H^/'^. 
Finally, the output modes are given by 0o"^(r)=0Q (x, ?/)0o^^^''(z), where n stands for the 
quantum numbers {kr^^ky, zeJ, and the density of modes is p^^e) = {1/87T'^)L'^H'^^'^ /I. 
Thus, the problem is reduced to the coupling of an initially populated bound state m=— 1, 



4 



of energy E-i=hSrf + /i to a quasi-continuum of final states m=0, with a total energy 
EQ^^=fi^{kl + ky)/2M — MgzE^ + 4yu/7. A crucial feature in this problem is the reso- 
nant bell-shape of the coupling matrix element W{i={fi^l^{/2){(f)Q^^ \ We work out 
the overlap integral /("^ = (0q"'^ | in reduced cylindrical coordinates {r~±,6,z) and set 
k"^ = {kxXoY + {kyUof' . We integrate over 9 and rj_, and transform the sum over z with the 
Parseval relation llSl to obtain: = 2nAxQyo/L{^i/U)^^/'^^ j^^g{k,v)e'^^^'^'^^'''' /^-"'^^/^Mv. 
In this integral, g is the Fourier transform with respect to z of g{k^ z) = {pcosp — smp)/k^, 
with p(k,z) = k{l — z'^y^'^, and (27r)'^"-'^/^)e*'^''^/^~^^^^/') the Fourier transform of Ai. Since 
zq/1 ^ 1, the integrand averages to zero out of a small neighbourhood of the origin, where 
the linear term in v is dominant. We obtain in this way: 

Al rji 



li-) = 27r-^ L^oyogCk,z~Ej (4) 

This overlap integral is non-vanishing only if the accessible final energies _i are 

restricted to an interval [—A/2, A/2], where A ~ 2MgzQ ~ 22.7 kHz. This gives a resonance 
condition for the frequency cUrf |P, ITSj : 

\h5rf\<Mgzo (5) 



Two different behaviour can be expected in such a situation [|T3], |T^. In the strong coupling 
regime {hV ^ A, where F is the decay rate worked out with the Fermi Golden Rule), Rabi 
oscillations occur between the BEG level and the narrow-band continuum. This describes 
the pulsed atom laser experimentally realized by Mewes et al. 0. On the contrary, in the 
weak coupling limit {hV ^ A), oscillations persist only for t < tc=h/A, while for t > tc, 
the evolution of the BEG level is a monotonous decay. We have numerically verified this 
behaviour on a ID simulation analogous to (see Fig.([l|a)). 

Eq.(|D can be formally integrated |13| with the help of the propagator Go of Hq: \I'o(r, t) = 
^ Jodt' J d^r' Go{r,t;r' ,t')'^ ^i{r\t'). Together with Eq.([l|) and Eq.(|^), we can derive an 
exact integro-differential equation on a^, the fraction of atoms remaining in the BEG. In 
the weak coupling regime, the condition /iF -C A expresses that the memory time of the 
continuum tc is much shorter than the damping time of the condensate level. This allows 
one to make a Born-markov approximation |TB|, which yields the rate equation: 

^ = -r(iv_i) iv_i (6) 
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The output rate T is given by the Fermi golden rule: 

with ri=2{h6r{ + 4/i/7)/A. The rate equation (H) is non linear, since T depends on A^_i 
through A = 2MgzQ. To obtain Eq.(^, we have neglected the transverse kinetic energy as 
compared to E^i, and used the relations fi/UxoyoZo = IS/Svr (valid in the TF approximation) 
and J(^°° dw \ wcos{w) — sin (w) p w"^=l/4. A quasi-continuous output corresponds to 
the weak coupling regime {hT <^ A). From Eq.(0), we deduce a critical Rabi frequency 
~ 0.8A/h for which hV ~ A. However, we have assumed from the beginning that 
the BEC decay was adiabatic. This requires | dh^i/dt |~ F/i ^ ^fi)/^^ where e(j) is the 
energy of the ith quasiparticle level in the trap. Taking e(j) > hu^, we deduce from Eq.(|^) a 
condition on the Rabi frequency: firf ^ 1.6{g/zQy^'^. This condition turns out to be much 
more stringent than the condition for weak coupling Qj.f ^ fi^. 

To compare our model to the data of Bloch et a/ 0| , we have numerically integrated Eq. 
with the output rate (0) and their experimental parameters. We show in Fig. (j^) the number 
of atoms remaining in the condensate after a fixed time as a function of the detuning. The 
agreement with the experimental data for the | F = l^mp = —1) state is good (Fig.(^)). 
To treat the | F = 2;mF = 2) case, the second trapping state | F = 2;mF = 1) must be 
included. We use the result of [jl9| for the situation where two discrete levels are mutually 
coupled, and only one of them is also coupled to a continuum (with a decay rate r2,i)[jl9|. 
For near resonant coupling, the upper level acquires a decay rate T2,2 ~ ^2,1 if ^rf ^ r2,i- 
We have taken this value as a first approximation. This leads to an already good agreement 
with the experimental data (Fig.^b), although better refinements is needed to describe the 
full dynamics. 

Finally, we want to point out that, since /iF ^ A, the spatial region where outcoupling takes 
place, of vertical extension 6z ~ hT/mg, is very thin compared to the BEC size. It allows us 
to think of the outcoupling process in a semi classical way, in analogy with a Franck-Condon 



principle [^: the coupling happens at the turning point of the classical trajectory of the 
free falling atoms. Neglecting the transverse kinetic energy in the calculation of F amounts 
to approximate the Franck-Condon surfaces by planes. Since over the size of the BEC, the 
surface curvature is small, the deviation from this plane is negligible compared to Zq. 
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We have analyzed so far the extraction of the atoms from the trapped BEC. We address now 
the question of the propagation of the outcoupled atoms under gravity. For times t ^ tc, 



the exphcit expression of the propagator Go of is 



Go{r,r',T = t-t) = (J^)3/2ei5,,,,,(r,r',r=t-t')g(^^ ^g) 

Here, the classical action is given by S'ciass(r, r', r = t — t') = (M/2/ir)((r — r')^ + gT'^{z + 
z) —g'^T^/ 12) +4/ir/7/i and 9 is the step function. We compute the output wavefunction \l/o 
with stationary phase approximations. This amounts to propagate the wavefunction along 
the classical trajectory: in the limit where S'dass ^ ^(this is true as soon as the atom has 
fallen from a height I), this is the Feynman path that essentially contributes. We obtain 
the following expression for the outcoupled atomic wave function, i.e. the atom laser mode 
(see Fig.(|l|b)): 

ipo{r,t) -^/7^—^(|)^l{x,y,Zres) F{t, tM\) (9) 

\ \ Cres 1^/2 



In this expression, Zj.es=i]Zo/2 is the point of extraction, tii^ii={2/ gY^'^^z — Zj-csY^"^ is the time 
of fall from this point, Cres=(-2— -Zres)/^ and F(/5)=a(/5)M(/5): M is the top hat function, that 
describes the finite extent of the atom laser due to the finite coupling time. We can deduce 
from Eq.@ the size of the laser beam in the x direction a;out = a^o(l — and 
a similar formula for yout- Because the semi-classical approximation neglects the quantum 
velocity spread due to the spatial confinement of the trapped BEC, this wavefunction has 
plane wave fronts. This property will not persist beyond a distance zji ^ y^^^/Al^ ~ a few 
mm, analog to the Rayleigh length in photonic laser beams, where the transverse diffraction 
{h/ Myont){2zR/ gy/"^ becomes comparable to the size ?/out- Beyond zr, diffraction has to 
be taken into account. However, this discussion neglects the (weak) transverse acceleration 
provided by the repulsive mean field potential that is likely to affect the transverse profile. 
In conclusion, we have obtained analytical expressions for the output rate and the output 
mode of a quasi-continuous atom laser based on rf outcoupling from a trapped BEC. Our 
treatment, which fully takes the gravity and the 3D geometry into account, leads to a good 
agreement with the experimental results of Ref . @] . This points out the crucial role played by 
gravity in the atom laser behaviour. More elaborated treatments should take the interaction 
energy between the trapped BEC and the atom laser, as well as finite temperature effects 
jlSl into account. Nevertheless, we believe that the present zero-temperature theory is a 



valuable starting point to describe experiments with the atom laser. 
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FIG. 1: Atom laser in the weak coupling regime (ilj-f = 300 Hz), la : time evolution of the laser 
intensity starting with ~ 2000 atoms. For t ^ tc, the numerical integration of Eq.(l) agrees with 
the output rate of Eq.(^). lb : Spatial intensity profile of the atom laser at t ^ tc according to 
Eq.d) 
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FIG. 2: Number of trapped atoms after 20 ms of rf outcoupling, starting with ~ 7.2 x 10^ condensed 
atoms, rirf=312 Hz for the | 1; -1) sublevel (left), and 7.0 x 10^ atoms, r2rf=700 Hz for the 
I 2; 2) sublevel (right). Diamonds are the experimental points from Bloch et al, solid line is the 
prediction based upon our model using the experimental parameters. Theoretical and experimental 
curves have been shifted in frequency to match each other, since VoS is not experimentally known 
precisely enough (a precision of 10~^ Gauss for the bias field Bq is required to know Vqs within a 
kHz uncertainty). 



10 




1 1 1 1 1 1~ T 1 1 1 1 r- 

-10 10 -10 10 



Detuning (kHz) Detuning (kHz) 



